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ABSTRACT 


UPDATEDSPOSITION AN DSICE VELOCITY 
FOR THE AIDJEX MANNED CAMPS 


A.S. Thorndike and T.O. Manley 


During 1975 and 1976, measurements of sea ice motion were made as 
part of the Arctic Ice Dynamics Joint Experiment (AIDJEX). In an earlier 
report (Thorndike and Cheung, 1977), energy at the inertial period was at- 
tenuated by approximately 30% due to the choice of a low covariance factor, 
Q = 100, during Kalman filtering. 

For oceanographic investigations, the data set proved to be unac- 
ceptable particularly in the summer months when there is a significant amount 
of energy at the inertial period. 

This report is the result of the reprocessing of the original position 
data taken from the four manned camps (Caribou, Blue Fox, Snowbird and 
Big Bear) with a covariance factor Q = 100000. 

The higher Q allows data of the inertial period and longer to pass 
through the filter with less than 1 percent attenuation. 

Because of the large amount.of information contained in the data set, 
it was divided into two parts. This report (Volume 1) provides position and 
ice velocity from 11 April 1975 to 17 October 1975. Volume 2 covers the re- 
mainder of the data, 18 October 1975 to 4 May 1976 (Thorndike and Manley, 


1980). 
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Introduction 


The objective of the Arctic Ice Dynamics Joint Experiment, as set 
forth in its most general statement, was to reach an understanding of the 
dynamic and thermodynamic interaction between the ice cover and its en- 
vironment (Maykut et al., 1972). One component of the experiment was an 
array of drifting stations at which frequent measurements of the position 
were made (Figure 1). For the present report those measurements have 
been edited using Kalman smoothing techniques to give the best estimates of 
position and velocity. (The raw measurements themselves are not reported 
here, but are available through the Polar Science Center.) The data set 
described here is restricted to the main AIDJEX experiment, which began in 
March 1975 and ended in May 1976. 

This data report represents the most complete documentation of ice 
motion available, but it is limited to one year and to the Beaufort Sea. Par- 
ties interested in earlier data, more widely dispersed in space and time, 
should contact the Polar Science Center, University of Washington, Seattle, 
Washington, which has compiled most of the ice drift data beginning with the 


drift of the FRAM in 1893 to the present. 
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Data Acquisition 


A full description of the measurement systems used in this program 
is beyond the scope of this report. The list of references includes sources 
of more detailed information (Wasilew and Vivian, 1976; Martin, P., 1978). 
All of the hardware we used was supplied by the Satellite Positioning Cor- 
poration, Houston, Texas. 

Four manned camps -- identified by radio call names Big Bear, Cari- 
bou, Blue Fox and Snowbird -- were tracked using NavSat (Navy Navigation 
Satelite System). The geographical position of each of the four camps was 
measured at unevenly spaced times using Doppler satellite navigation tech- 
niques. The raw measurement is the frequency of a signal transmitted from 
the satellite to the ice station. The measured frequency is affected by a 
Doppler shift related to the rate of change in distance between the satellite 
and the ice station, which is itself related to the unknown ice coordinates 
and the known satellite orbit. During each satellite orbit several frequency 
measurements were made and from them an overdetermined solution was found 
for the unknown coordinates. Data recording in the field was on magnetic 
tape with hard copy printout for backup. 

A constelation of six polar orbiting satellites was maintained for this 
system. The orbital period is about 107 minutes. An observer at high lati- 
tudes can track every orbit. Of the 81 potential passes each day, many are 
unusable (because of radio interference between overlapping passes, for in- 
stance). We averaged about 25-30 usable passes per day at each station. 


Single channel (400 MHz) receivers were used to count the Doppler cycles 


for roughly 30-second intervals. The duration of each Doppler counting 
interval was measured independently. The frequency standard used in the 
receivers was stable during a 20-minute pass to about 1 part in ne 

Fix solutions were done in real time at each manned camp. The 
solution used only the Doppler data from the 5 minutes on either side of the 
time of closest approach of the satellite. No additional editing was performed. 
A correction to each Doppler count was based on the duration of the counting 
interval to account for the random variations in the counting interval. The 
correction also included the fixed timing delay in each receiver. The fix 
solution ignored the effect of the velocity of the ice camps. 

In fact, two receivers with antennas separated by about 100 meters 
were operated at each manned camp. In addition to serving as a backup, 
the measurements at the second antenna were used to define the azimuth of 
the line joining the two antennas. This was possible when both receivers 
produced acceptable fixes from the same satellite pass. Estimates of the 
azimuths obtained at each of the manned camps were reported by Thorndike 


and Cheung 71977. 


Preprocessing 


The original data tapes from the manned stations measurement sys- 
tems typically contained several "glitches" per tape because of power failures 
at the camp, or operator error, or some other reason. These tapes were 
transcribed discarding those records which had either a parity error or an 
incorrect number of words. The tapes which resulted from this process 
were read without error on the CDC 6400 at the University of Washington, 
where the bulk of the processing was done. 

Three problems had to be dealt with to prepare the sequence of 
measurements at each station for the actual Kalman filtering. First, data 
recording problems had prentad occasional gaps in the manned camp data. 
At times the backup teleprinter records from the field could be used to fill 
these gaps. This was done by keypunching the fixes directly from the 
printed records. Second, the data sequence had occasional points which 
appeared to be out of sequence in time. This was sometimes correctable, 
as when the operator had keyed in the wrong date, but usually the situation 
was ambiguous and points out of chronological order were dropped. The 
third and most serious problem involved bad fixes. A glance at a plot of 
the raw time series (Figure 2) illustrates the problem. Because outlying 
fixes have a deleterious effect on the Kalman filtering results, an attempt 
was made to remove them during preprocessing. The algorithm employed 
compares each fix in latitude with the median of the latitude of the 10 fixes 
preceding and the 10 fixes following it. When the differences exceeded a 


preset tolerance of 10 kilometers, the fix was eliminated. The algorithm 
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was applied twice, with a smaller tolerance of 1 kilometer the second time. 
The same procedure was then applied for longitude. 

Visual checks were an important part of the procedure. Plots of 
the data points were produced before and after preprocessing. At times it 
seemed from the plots that we were not detecting certain bad fixes, or that 
we were throwing out good ones, and it was necessary to adjust the toler- 
ance levels. Edited raw data plots like the one shown in Figure 3 provided 
assurance that the data were ready for filtering. 

For the main experiment, AIDJEX adopted the convention of number- 
ing days consecutively, beginning with day 1 = 1 January 1975 and ending 
with day 500 = 14 May 1976. Throughout this data waREe, this numbering 
system is used quite extensively. Appendix 3 provides a conversion table 


from AIDJEX days to calendar days. 
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Kalman Smoothing 


In the smoothing used here an assumption is made regarding the 
motion of the ice. If the state of the ice at time th is represented by a 
vector Xn? containing, say, the position, velocity, and acceleration of a 
point, then the state at some future time is assumed to be partially deter- 


mined by ae 


Kevee= OX el Wo (1) 


where © and T are matricies known from the ohysics of the situation, and 


WwW represents an unknown random perturbation. In Kalman's formulation, 


n+l 


which we follow, the random effects are assumed to be Gaussian and White: 


Q ifn=m, and 
cov (W_W_) = (2) 
ha 0ifn¢m 
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wes 
We have chosen Q = 100000 m°hr ~. Qualitatively choosing a small value of 


Q is equivalent to assuming that the ice moves in a completely predetermined 
way. (See from equations 1 and 2 that Q = 0 would imply that A ntl = 9 xn 
a large value of.Q corresponds to assuming that the ice experiences large 
sudden changes in its acceleration.) Thorndike and Cheung (1977) used 

Q = 100 mona The lower value of Q assumed smooth behavior, but was 
unnecessarily conservative. The present value of Q was determined by trial 
and error. aA sequence of raw fixes was used as input to the Kalman filter. 
Several different runs were done using different values of Q. The present 


value of Q appeared to give the best agreement between the output from the 


filter and the raw fixes. 


The original Q (100 m* hr ®) used by Thorndike and Cheung removed 
approximately 50% of the inertial period amplitude as can be seen in Figure 4. 
A constant amplitude sinusoidal function whose period varied linearly with 
time was chosen as the input to the Kalman filter with Q = 100 chr eee 
put periods ranged from 0.5 hours (at day 10) to 60.4 hours (at day 15) with 
the inertial period being located at day 10.9. An approximate 50% attenuation 
of the input function at the inertial period can be roughly estimated from the 
graph. 

Figure 5 shows the same input function and the resulting Kalman 


28, The inertial period passed through 


output with Q chosen to be 100000 m 
the filter with very little attenuation. Another test was then run, this time a 
constant amplitude, constant period (12 hours) sine wave was used to simulate 


the inertial period of the Arctic latitudes. The results indicated that there 


was only a 0.03% attenuation from the original input function. 


For the Kalman reduction, we defined A = te = ta? and the 
matrices 
pap ; 2 
x = position d = Va saat 2 
velocity De ae (3) 
acceleration 7550 0 
ree “aKa 
S 
Net2 (4) 
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ihe dimensions of the matrices are X :3x1,9Q:3x3, [:3x1, 0:3x3, 
Wee cali. erate . 

An identical equation can be written for the y-coordinates. The 
two coordinates are treated independently. 


The measurement Z4 is related to X by 
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; ; é 2 
where Meee is a Measurement error with assumed covariance R = (60 meters) . 


The response of the measurement and processing system was investi- 
gated as follows. Four hundred sample points were drawn independently 
from a Gaussian distribution with zero mean and variance (50 meters) ?. A 
sampling interval of one hour was assumed. This noise simulated the Nav- 
Sat measurements. The sequence was used as input to the Kalman filter 
procedure. The velocity output from the filter shows two effects. First, 
since velocity is the time derivative of position, we expect to find the spec- 
trum S,,(e) behaving as w"S (uw); where S, Cw) is the spectrum of position. 


2 yey. ey eee 
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S$, (4) was assumed to be white; S,. Cw) 
The second effect is the tendency of the filter to smooth out high frequen- 
cies. Qualitatively then, we have something like the diagram in Figure 6. 
The actual spectrum of ice velocity is shown in Figure 7. The dia- 
gram illustrates how the measurement errors and finite sampling rate limit 
the accuracy of the velocity estimates. Most of the variance in the velocity 
errors occurs for the frequencies greater than 3 cycles per day. For lower 


frequencies we can estimate the contributions to the variance as 0 - 1.3 cycles 
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Figure 6 - Diagram showing qualitatively, the response of 
the measurement and processing system. 
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processing system. For a more complete de- 
scription, see text. 
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per day: o7 = 3x10 °m “le =, 1.5- 3.0 cycles per day: o7 =2x10° 
mn sec - total variance 0 - 12 cycles per day: a” =e ghe m* sec. 

The general smoothing problem is to estimate X at time Me using the 
entire set of observations (Z;) ijsul,.-.M, and to give the error variance 
of that estimate. The solution equations, due to Kalman, are reproduced 
in Appendix 1. All positions were converted to an X-Y coordinate system 
prior to actual filtering. This coordinate system along with its governing 
equations are shown in Appendix 2. Typical estimation error variances for 


position and velocity are 40 meters and 1.0 cm/sec respectively. 


Results 


The final data listing consists of five major columns. The first 
column is the time in decimal AIDJEX days corresponding to the observation 
of the position and ice velocity at a particular camp. The remaining columns 
indicate the camps for which data are applicable. Subdivisions under each 
major camp heading are LAT (latitude in decimal degrees), LONG (longitude 
in decimal degrees), and the north and east components of ice velocity in 
em/sec, NIV and EIV respectively. 

Due to the large amount of information contained in the data set, it 
was divided into two parts. This report (Volume 1) covers AIDJEX days 
101 (11 April 1975) to 290 (17 October 1975). Volume 2 covers AIDJEX days 
291 (18 October 1975) to day 490 (4 May 1976) (Thorndike and Manley, 1980). 

Figures 8 - 11 show a more detailed drift track of each of the manned 
camps. An insert shows where the region lies with respect to the north slope 
of Alaska and Banks Island. Asterisks superimposed on the drift track indi- 
cate the positions at integral multiples of 20 days. The beginning and end- 
ing days are noted for each trajectory. A dashed line indicates a long period 
of missing data. The plotted region is 500 x 500 km and is aligned with the x, 
y coordinate system shown in Appendix 2. Figures 8 - 1l were reproduced 


from Thorndike and Cheung (1977). 


16 


MoSEI 


M,OvT 


MoSYT 


Nov 


G9) 


N94 


N,SZ 


Notl 


MSI 


SII 


“28b USsnoOuy) 
Aep woay 1RG AN 


Aq sjyuswa.inseoay, «dn 


ayo 


iq avog Stgq aye 


dweod urew ay) ouweoog 


MOST 


noqmeg dwuer9g 


~ § dIndty 


17 


MNoOLT 


MoSET 


M071 


OOS 


MoS7I 


NoSL 


Maal 


“6LG 

Aep uo woay p.reoy sep sem 
Aong Uy, “€LG 01 LEG Woay 
ded vyep 8 St ota, “posn 
DIOM ZLLUW AON SINVUL WoOuy 
e1ep 06Gb Arp aay “198s 
B1ep oy UT pasn jou 9.19M 
nq potaod yey. Aurinp spew 
aIOM SOXTT uotpisod [rnsojag 
"802 01 €0z Woay ded eB uM 
68h YUdnoOAYI OTT ARep woay 
WS AGN Aq SUDWOANSR|N 


xO ont” dweg 


~ § dandy 


MOGI 


NS Niet N,€Z MaGill 


x Z MoS 
Do Sicat 
*SOUWT] [RUDADS 
N.OvL SOOURISIP 1.LOYUS PAO 210M 
SBUUOJUP TEGAN 9Y1 pur 
IdM OY] JULINp syora9 
Aq 1uds sem uorneis sry, 
“T6b UdNOAYT GOT wo 
wsaen Aq sjuoWwDuNsrasy 
parqmoug dung 
~ OL aan 
SISAL 
MoSGT 
NLOSL K 


NoLL No9L NGL NoVL 


ao 


Nov No€L Nocl 


“SINTOWOTHy pPRa19 
-~AdS TDAO peroaids u90q 
pey dwro pjo oy} own 
sty1 Ag “Sgp-sop sAup 
‘poou Aong SWVH dutsn 
apew 919M SPUDWOANS POW 
uortsod pue pordnsso0 
-d91 SPM UOTIRIS 9U1 OLGT 
jo duridg oy ur potaad 
jJoriq & 107 "087 ABD 

uo dulaq eyep yseR_ oud ca 
‘pauopurqe ATpenjUudAd 
pue Aytarjor vor Aq 
dn usyo1qg sem dwro ay 
bLeZ Avp uo “201 Avp 
uo duruutdaq yegaen Aq 
APRW IAIM SPUIWIANS IIA] 


anog dig dung 


Tt oansdry 


MSI 


Noll No9/ NoSL 


APPENDIX 1 


KALMAN SMOOTHING EQUATIONS 


Symbols 

pee Cae Vy eClLOrat x11) 

@ = transition matrix (n xn) 

[ = coefficient matrix for random perturbations to physical 
model (nx m) 

W = white noise driving the physical process (mx 1) 

Q = covariance matrix for W (m xm) 

Z = measurement vector (2 x 1) 

H = # matrix relating measurement to state vector (2 xn) 

V = measurement error vector (2x 1) 


= measurement error covariance (2 x 2) 


P = prediction, filtering and smoothing error covariance 
matrix (n x n) 

I = identity matrix (n xn) 

K = £Kalman gain matrix for filtering (n x 2) 

A = smoothing gain matrix (n x n) 

N = number of measurements in the data sequence 

kK = index for the data sequence 

c = time 


The superscript T denotes the transpose operation. The symbol 
X (i/j) means the estimate of y at time t given observations (Z,, Zo; ane 7 aa 
The notation here follows Meditch (1969). 


Prediction 


Given observations (Z,, as Z,), predict the state at the future 


time t, and find the prediction error covariance. 


Koen 


wick al, | kK) re (Keele eee, A(1) 


P(k+1/k) = 6(k+1,k)P(k/k) 6? (k+1,k) 
+ Tk +1,k)Q(k) Pr? (k= 1,k). A(2) 
Filtering 
Given observations (2, sree Z,.) and a new observation Z +1)? 
estimate the present state at time ted and find the filtering error co- 
variance. 
X(k+t1[k+1l) = Lck+1[k) +K (K+ D(ZK+ 1) 
- H(k +1) xX(k+1/k)), A(3) 
P(ko} 1) +1)? = eK ee ee a eee A(4) 
where 


K(k +1 = P(k+1[k) Bw (k+ D (H(K +1) P(e + 1/kK) BP (Ck + D 


Reve Yate a aahie A(3) 


Smoothing 


Bie et eth 


Given the complete data set (Z,, i sels 2) , estimate the state at 


any point thee and find the smoothing error covariance. 


x (k | N) 


Y(k|k) +A (kK) (X(K +1])N) - X (kK +1/k)) A(6) 


P(K|N) = P(k|k+A(k) (P(K+1|N) 
- P(k+1/k)) A! (k), where A(1) 


A(k) = P(k|k) 67 (k +1,k)P 1 (Kk + 1/k). A (8) 


to 
bo 


Implementation of the Filter Algorithm 

Starting from an initial condition x (0/0), P(0| 0), a prediction is 
made to time ty. The prediction error covariance matrix P (1! 0) is used to 
obtain the Kalman gain K(1). The filtered estimate ~(1| 1) is then found 
by adding to the predicted yx (1| 0) a quantity which is proportional to the 
difference between the observation Z(1) and H(1)¥(1/0). H(1)%(1]0) is 
the value the measurement would have if the true state were actually %(1/ 0) 
and no measurement errors were made. The algorithm is repeated as each 


new data point is entered until the last filtered quantities % (N|N) and 


P(N |N) are produced. 

The smoothing algorithm, proceeding backward in time, utilizes 
results from the prediction and filtering calculations. The first step is to 
find { (N-1|N) and P(N-1|N) from equations A(6), A(7) and A(8). No 
direct use of the measurements is made in smoothing. The algorithm simply 
corrects the filtered vector xy (kK|k) by a quantity proportional to the differ- 
ence between the prediction vector at time tad and the smoothed value at 
time teed: 

In the calculations done for this application, the CDC 6400 computer 
time required to do the prediction, filtering, and smoothing is about 0.3 sec- 
ond per data point, nearly half of which is absorbed in the matrix inversion 


in A(8). Because of the recursive nature of the filter, demands on central 


memory storage space are very modest. 


Reproduced from Thorndike, A.S., 1973, An integrated system for meas- 
uring sea ice motions. In Ocean 73, IEEE International Conference on En- 
gineering in the Ocean Environments, IEEE publ. no. 73 CHO 774-0 OEC, 
Dore I0-490., 
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APPENDIX 2 


COORDINATE SYSTEMS 


Position measurements were made in geographical coordinates 
(latitude north, longitude east). The smoothing operation was 


done in a Cartesian system (x,y), where 


x = 110.949 (90° - latitude) cos (longitude) (km) 


as 
i 


110.949 (90° - latitude) sin (longitude) (km) 


APPENDIX 3 


CONVERSION TABLE FOR AIDJEX DAYS 
TO CALENDAR DAYS 


For the main experiment, AIDJEX adopted a convention of number- 
ing days consecutively, beginning with day 1 = 01 January, 1975 and ending 
with day 500 = 14 May, 1976. 

In the conversion table, the first column is the AIDJEX day, the 
second is the corresponding day of 1975 or 1976 and the third entry is the 


ealendar date. 


Appendix 3 - Conversion Table 
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Position and Ice Velocity Data 


The following section comprises the updated positions and ice veloci- 
ties for the four manned camps Caribou, Blue Fox, Snowbird and Big Bear. 


Data span from 11 April 1975 to 17 October 1975. 
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